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The formation of a naked singularity in a model of f(R) gravity having as source a linear elec- 
tromagnetic field is considered in view of quantum mechanics. Quantum test fields obeying the 
Klein— Gordon, Dirac and Maxwell equations are used to probe the classical timelike naked sin- 
gularity developed at r = 0. We prove that the spatial derivative operator of the fields fails to 
be essentially self-adjoint. As a result, the classical timelike naked singularity remains quantum 
mechanically singular when it is probed with quantum fields having different spin structures. 
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I. INTRODUCTION 



In the last decade, there have been extensive studies in Extended Theories of Gravity (ETG) such as the Lovelock 
and f(R) gravity theories. The main motivation to study the ETG is to understand the accelerated expansion of 
the universe and the issue of dark matter/energy (see [l[ and references therein for a general review). One of the 
I/*-) ' most attractive branches of the ETG is the f(R) gravity theory in which the standard Einstein's gravity is extended 
with an arbitrary function of the Ricci scalar R instead of the linear one In this model, the Ricci scalar R in 
the Einstein— Hilbert action is replaced with f(R) = R + ag(R), where g(R) is an arbitrary function of R so that in 
' ^j' , the limit a = 0, one recovers the Einstein limit. Although the majority of researchers prefer to use this ansatz, in 
qh> general, finding an exact analytic solution to the field equations is not an easy task. As far as analytic exact solutions 
are concerned, static, spherically symmetric models in f(R) gravity have been shown to serve for this purpose 041] • 
(3J[)' In this context of static, spherically symmetric solutions of f(R) gravity, the solutions admitting black holes have 
attracted much attention. 

In the context of static, spherically symmetric f(R) gravity, it has recently been shown that 0], an exact analytic 
solution is also possible if one assumes f(R) to have the form of f(R) = £ (R + Ri) + 1ol\JR + Rq, in which £, a, Rq 



ly-j , and i?i are constants, a priority to secure the Einstein limit by setting the constants Rq = Ri = a — and £ = 1. 

■ In this model of f(R) gravity, exact solutions with external electromagnetic sources (both linear and nonlinear) are 
7— I | found. It was shown that the solution with a linear electromagnetic field does not admit a black hole while the solution 
with a nonlinear electromagnetic source admits a black hole solution. The physical properties of the latter solution 
| are investigated by calculating thermodynamic quantities and it was shown to satisfy the first law of thermodynamics. 
. The solution having as a source a linear electromagnetic field resulted with a naked curvature singularity at r = 0, 
C^l 1 which is a typical central singularity peculiar to spherically symmetric systems. The solution given in Q, is a kind 
of extension of a global monopole solution Q which represents a solution of the Einstein's equations with spherical 
1 symmetry with matter that extends to infinity. It can also be interpreted as a cloud of cosmic strings with spherical 
symmetry Q. Hence, the spacetime is conical. However, with the inclusion of a linear or nonlinear electromagnetic 
rS ■ field, the spacetime is no more conical in the context of f(R) gravity. 

Within the framework of ETG gravity, black hole solutions have been widely studied in the literature (see 0, Q3 and 
references therein for a complete review). However, the solutions that result with naked singularities have not been 
studied in detail. In physics, naked singularities are considered to be a threat to the cosmic censorship hypothesis. 
Furthermore, as in classical general relativity, compared to the black hole solutions, naked singularities are not well 
understood in the context of f(R) gravity. This still remains a fundamental problem in general relativity as well 
as in ETG to be solved. Another important difficulty in resolving this problem is the scale on which the curvature 
singularity occurs. On these small scales, it is believed that the classical methods should be replaced with quantum 
techniques in resolving the singularity problems that necessitate the use of quantum gravity. Since the quantum 
theory of gravity is still "under construction", an alternative method is proposed by Wald [TiJ which was further 
developed by Horowitz and Marolf (HM) [12| in determining the character of classically singular spacetime and to see 
if quantum effects have any chance to heal or regularize the dynamics and restore the predictability if the singularity 
is probed with quantum particles/fields. 
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In this paper, we investigate the occurrence of naked singularities in the context of f(R) gravity from the point 
of view of quantum mechanics. We believe that this will be the unique example wherein the formation of a clas- 
sically naked curvature singularities in f(R) gravity will be probed with quantum fields/particles that obey the 
Klein— Gordon, Dirac and Maxwell equations. The criterion proposed by HM will be used in this study to investigate 
the occurrence of naked singularities. 

This criterion has been used successfully for other spacetimes to check whether the classically singular spacetimes 
are quantum mechanically regular or not. As an example; negative mass Schwarzschild spacetime, charged dilatonic 
black hole spacetime and fundamental string spacetimes are considered in Il2l| . An alternative function space, namely 
the Sobelov space instead of the Hilbert space, has been introduced in [13[ , for analyzing the singularities within the 
framework of quantum mechanics. Helliwell and Konkowski have studied quasiregular [lj], Gal'tsov— Letelier— Tod 



spacetime [151 ] , Levi-Civita spacetimes [ly, [17[ , and recently, they have also considered conformally static spacetimes 
[lq . Pitelli and Letelier have studied spherical and cylindrical topological defects Banados— Teitelboim— Zanelli 
(BTZ) spacetimes flcj, the global monopole spacetime [2l[ and cosmological spacetimes [22j . Quantum singularities 
in matter coupled 2 + 1 dimensional black hole spacetimes are considered in |23| . Quantum singularities are also 
considered in Lovelock theory [13] and linear dilaton black hole spacetimes [IJ Recently, the occurrence of naked 
singularities in a 2 + 1 dimensional magnetically charged solution in Einstein— Power— Maxwell theory have also been 
considered (26|. 

The main theme in these studies is to understand whether these classically singular spacetimes turn out to be 
quantum mechanically regular if they are probed with quantum fields rather than classical particles. 

The solution to be investigated in this paper is a kind of f(R) gravity extension of the analysis presented in [2l[ 
for the global monopole spacetime. The inclusion of the linear Maxwell field within the context of f(R) gravity 
affects the topology significantly and removes the conical nature at infinity. Furthermore, the true timelike naked 
curvature singularity is created at r = which is peculiar to spherically symmetric systems. We investigate this 
singularity within the framework of quantum mechanics by employing three different quantum fields/particles obeying 
the Klein— Gordon, Dirac and Maxwell fields with different spin structures. 

The paper is organized as follows: In Sec. II, we review the solution found recently in Q, and give the structure of 
the spacetime. In Sec. Ill, first, the definition of quantum singularity for static spacetimes is briefly introduced. Then, 
the quantum fields obeying the Klein— Gordon, Dirac and Maxwell equations are used to probe the singularity. The 
paper ends with a conclusion in Sec. IV. 

II. THE METRIC FOR f{R) GRAVITY COUPLED TO MAXWELL FIELDS AND SPACETIME 

STRUCTURE 

Recently, an exact analytic solution for / (i?) gravity coupled with linear and nonlinear Maxwell field in four 
dimensions has been presented in [tJ ■ The corresponding action for / (R) gravity coupled with linear Maxwell field in 
four dimensions is given by, 



fM_± F 

2k 4tt 



(1) 



in which / (R) is a real function of the Ricci scalar R, and F — \F^ LV F^ V is the Maxwell invariant. The Maxwell 
two-form is given by 

F = Qdt A dr + P sin 9d9 A dip, (2) 

in which Q and P are the electric and magnetic charges, respectively. The static spherically symmetric metric ansatz 
is 

dr 2 

ds 2 = -B (r) dt 2 + — — — - + r 2 (d9 2 + sin 2 Odip 2 ) , (3) 
B yr) 

where B (r) stands for the only metric function to be found. The Maxwell equations (i.e. dF = = d*F) are satisfied 
and the field equations are given by 

j r r; + (nf R - M s; - v"v,/ fl = «t*, (4) 
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in which 



fR = ~diT' (5) 

□/fl = -^^(^"j/K, (6) 



-0 

v F vj« = <r [(/ B ), Uia - r™ (/ H ) tm , (?) 



while the energy momentum tensor is 



Furthermore, the trace of the field equation (4) reads 



AnTY = -FF ll + F ll xF uX . (8) 



2 

with T = Tjf. The non-zero energy momentum tensor components are 



f R R+(d-l)nf R -~f = K T, (9) 



M 87rr 4 

and with zero trace, we have 



= ~~Z~~T~diag [—1, — 1, 1, 1] , (10) 



2 

With reference to the paper [7j , the form of the function / (R) is assumed to be , 



f = lhR + -iaf R . (ii) 



/ (R) = £ (r + ii? ) + 2aV^ + Ro, (12) 



which leads to 



R — 2 r 2 ~ 

where a , Ro, and £ are constants. Consequently, the metric function B(r) is obtained for the free parameters a — rj 
as, 

„ , , 1 m o 2 A P ff „ , . 

B(r) = --- + q --^r 2 , 14 
2 r r z 3 

where m = A e // = —j 2 - and g 2 = - ^ p . As was explained in Q, due to the constraints on the free parameters, 
this solution does not admit the Reissner— Nordstrom (RN)— de Sitter (dS) limit. However, in the limit £ = 1 and 
P = Q = 0, the solution reduces to the well known global monopole solution reported in which represents a 
spherically symmetric, non- asymptotically flat solution with a matter field that extends to infinity. Furthermore, this 
solution can also be considered as a spherically symmetric cloud of cosmic string which gives rise to a deficit angle 
Q. Therefore, the solution given in equation (14) , is a kind of Einstein— Maxwell extension of the global monopole 
solution in / (R) gravity. One of the striking effects of the additional fields is the removal of the conical geometry 
of the global monopole spacetime. The Kretschmann scalar which indicates the formation of curvature singularity is 
given by 

1 8A 2 r 8 + 4Ar 6 + 3r 4 + I2mr 3 + 12r 2 (3m 2 - q 2 ) - U4mq 2 r + +168<? 4 

^=3 7 s • 

It is obvious that r = is a typical central curvature singularity. This is a timelike naked singularity because the 
behavior of the new radial coordinate defined by r* = J -g^pj is finite when r — > 0. Hence, the new solution obtained 
in [3] and given in equation (14) is classically a singular spacetime. 

Our aim in the next section is to investigate this classically singular spacetime with regard to the quantum me- 
chanical point of view. 



4 



III. QUANTUM SINGULARITIES 

One of the important predictions of the Einstein's theory of general relativity is the formation of spacetime sin- 
gularities. In classical general relativity, singularities are defined as the points in which the evolution of timelike or 
null geodesies is not defined after a proper time. According to the classification of the classical singularities devised 
by Ellis and Schmidt , scalar curvature singularities are the strongest ones in the sense that the spacetime cannot 
be extended and all physical quantities, such as the gravitational field, energy density and tidal forces, diverge at 
the singular point. In black hole spacetimes, the location of the curvature singularity is at r = and is covered by 
horizon(s). As long as the singularities are hidden by horizon(s), they do not constitute a threat to the Penrose cosmic 
censorship hypothesis. However, there are some cases that the singularity is not hidden and hence, it is naked. In 
the case of naked singularities, further care is required because they violate the cosmic censorship hypothesis. The 
resolution of the naked singularities stands as one of the most drastic problems in general relativity to be solved. 

Naked singularities that occur at r — are on the very small scales where classical general relativity is expected to 
be replaced by quantum theory of gravity. In this paper, the occurrence of naked singularities in f(R) gravity will 
be analyzed through a quantum mechanical point of view. In probing the singularity, quantum test particles/fields 
obeying the Klein— Gordon, Dirac and Maxwell equations are used. In other words, the singularity will be probed 
with spin 0, spin 1/2 and spin 1 fields. The reason for using three different types of field is to clarify whether or not 
the classical singularity is sensitive to the spin of the fields. 

Our analysis will be based on the pioneering work of Wald, which was further developed by HM to probe the 
classical singularities with quantum test particles obeying the Klein— Gordon equation in static spacetimes having 
timelike singularities. According to HM, the singular character of the spacetime is defined as the ambiguity in the 
evolution of the wave functions. That is to say, the singular character is determined in terms of the ambiguity when 
attempting to find a self-adjoint extension of the operator to the entire Hilbert space. If the extension is unique, it is 
said that the space is quantum mechanically regular. A brief review now follows: 

Consider a static spacetime (M,g^) with a timclike Killing vector field Let t denote the Killing parameter and 
£ denote a static slice. The Klein— Gordon equation in this space is 



VV,, - M 2 ) ip = 0. (15) 



This equation can be written in the form 



o2 / 

■qF = y/fD* (V?AV) - fM 2 Tp = -Aip, (16) 

in which / = — and Di is the spatial covariant derivative on £. The Hilbert space H, (L 2 (£)) is the space of 
square integrable functions on £. The domain of an operator A, D(A), is taken in such a way that it does not enclose 
the spacetime singularities. An appropriate set is C§° (£), the set of smooth functions with compact support on E. 
The operator A is real, positive and symmetric; therefore, its self-adjoint extensions always exist. If it has a unique 
extension Ae, then A is called essentially self- adjoint |27H29l |. Accordingly, the Klein— Gordon equation for a free 
particle satisfies 



(17) 



with the solution 



ip(t) = exp -ity/As V(0). (18) 

If A is not essentially self-adjoint, the future time evolution of the wave function (18) is ambiguous. Then the HM 
criterion defines the spacetime as quantum mechanically singular. However, if there is only a single self-adjoint 
extension, the operator A is said to be essentially self-adjoint and the quantum evolution described by Eq.(18) is 
uniquely determined by the initial conditions. According to the HM criterion, this spacetime is said to be quantum 
mechanically non-singular. In order to determine the number of self- adj oint extensions, the concept of deficiency 
indices is used. The deficiency subspaces N± are defined by ( see Ref. [13[ for a detailed mathematical background) 



N + = {ip € D(A*), A*ip = Z+ifj, ImZ + > 0} with dimension n+ 
N- = {ijj e D(A*), A*i; = Z-ip, ImZ- < 0} with dimension n_ 



(19) 
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The dimensions ( n+, n_) are the deficiency indices of the operator A The indices n + (n_) are completely independent 
of the choice of Z + (Z_) depending only on whether or not Z lies in the upper (lower) half complex plane. Generally 
one takes Z + = iX and Z_ = — iX , where A is an arbitrary positive constant necessary for dimensional reasons. The 
determination of deficiency indices is then reduced to counting the number of solutions of A*ip = Zip ; (for A = 1), 

A*ip±iip = (20) 

that belong to the Hilbert space T-L. If there are no square integrable solutions ( i.e. n + = n_ = 0), the operator 
A possesses a unique self-adjoint extension and essentially self-adjoint. Consequently, the way to find a sufficient 
condition for the operator A to be essentially self-adjoint is to investigate the solutions satisfying Eq. (20) that do 
not belong to the Hilbert space. 

A. Klein Gordon Fields 

The Klein— Gordon equation for a scalar particle with mass M is given by 

□V> = g- 1/2 d^ [(^VU] 1> = M 2 ip. (21) 
For the metric (3), the Klein— Gordon equation becomes 

d 2 ip n/ J n/ ,^ 1 d 2 ip 1 d 2 ip cot 9 dip (2B(r) / . A dip 

+B (r) M 2 ip. 

In analogy with equation (16), the spatial operator A for the massless case is 

, n , J n „ 9 2 Id 2 Id 2 cot 8 d (2B(r) „, , A d 1 

and the equation to be solved is (^4* ± i) ijj = 0. Using separation of variables, rp = R(r) Y™ (6, ip), we get the radial 
portion of equation (20) as 

* R U , ^B(r))' dRi r) , HP + 1) ± M a(r)=a (24) 



(22) 



dr 2 r 2 B{r) dr \ r 2 B{r) B 2 ( 

where a prime denotes the derivative with respect to r. 

1. The case of i — > oo 

The case r — > oo is topologically different compared to the analysis reported in [21] . In the present problem the 
geometry is not conical. The approximate metric when r — > oo is 

ds 2 ~ ~{ R ^-)dt 2 + ( dr 2 + r 2 U6 2 + sin 2 9d<p 2 ) . (25) 



0, (26) 



12 1 \R r , 

For the above metric, the radial equation (24) becomes, 

d 2 R{r) AdR(r) 
dr 2 r dr 
whose solution is 

where C\ and C2 are arbitrary integration constants. It is clearly observed that the above solution is square integrable 
as r — > 00 if and only if C\ = 0. Hence, the asymptotic behavior of R{f) is given by R{r) ~ 
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2. The case of r-> 



Near the origin there is a true timelike curvature singularity resulting from the existence of charge. Therefore, the 
approximate metric near the origin is given by 

ds 2 ~ -(^)dt 2 + (^j dr 2 + r 2 {dO 2 + sin 2 9d^ 2 ) . (27) 

The radial equation (24) for the above metric reduces to 



dr 2 q 2 
whose solution is 



^ Z(Z + 1 V)=0, (28) 



R(r) = C 3 e ar + C 4 e- ar (29) 

Vi (i + i) 

a = — 

q 

where C3 and C4 are arbitrary integration constants. The square integrability of the above solution is checked by 
calculating the squared norm of the above solution in which the function space on each t = constant hypersurface £ 
is defined as % —{R \ \\ R\\< 00}. The squared norm for the metric (27) is given by 



II Rf= f" '*<f ,30) 
Jo 9 

Our calculation has revealed that the solution above is always square integrable near r = 0, even if I = 0, which 
corresponds to the 5-wave solutions. 

Consequently, the spatial operator A has deficiency indices n + = n_ = 1, and it is not essentially sclf-adjoint. 
Hence, the classical singularity at r — remains quantum mechanically singular when probed with fields obeying the 
Klein— Gordon equation. 



B. Maxwell fields 



The Newman— Penrose formalism will be used to find the source-free Maxwell fields propagating in the space of 
f(R) gravity. Let us note that the signature of the metric (3) is changed to —2 in order to use the source-free Maxwell 
equations in Newman— Penrose formalism. Thus, the metric (3) is rewritten as, 

dr 2 

ds 2 = B (r) dt 2 - - r 2 (d6 2 + sin 2 6d<p 2 ) . (31) 

B (r) 

The four coupled source-free Maxwell equations for electromagnetic fields in the Newman— Penrose formalism is given 

by 

Dfa-tyo = {n-2a)(j> + 2p ( j> 1 -K<i> 2 , (32) 

5<j) 2 -Afa = -Wo + 2W>i + ("T -2/3) 2 , 

= (n - 27) <p + 2rcj) l - cr02, 

D^-'S^ = -A0 O + 271-^ + (p - 2e) fo, 

where B(r) is the metric function given in Eq.(14), (fi , <p 1 and <f> 2 are the Maxwell spinors, e, p, tt, a, ji, 7, /3 and t are 
the spin coefficients to be found and the bar denotes complex conjugation. The null tetrad vectors for the metric (31) 
are defined by 

ia - (wv w )> (33) 

i -^0 

2' 2 



mr 



y/2 \ r rsmv J 
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The directional derivatives in the Maxwell's equations are defined by D = l a d a , A = n a d a and S = m a d a . We define 
operators in the following way 

D = D, 

D = -5^A, (34) 



B(r) 

Lj = V2r S and h\ = L.J 
Lq = V2r 5 and Li = Lq 



cot6> 
cot6> 



The non-zero spin coefficients are 



lB(r) 1 , 1 cot 6» 



The Maxwell spinors are defined by [3C 



o = F 13 = F^m" (36) 
f 

F 42 = F^WW, 



H = ^(Fi2 + F 43 ) = -^(IV+m^m"), 



where Fij (i, j = f , 2, 3, 4) and F^,, (/x, i/ = 0, 1, 2, 3) are the components of the Maxwell tensor in the tetrad and tensor 
bases, respectively. Substituting Eq.(34) into the Maxwell's equations together with non-zero spin coefficients, the 
Maxwell equations become 

Do + 0i - ^ L ^o = 0, (37) 
Do + ^ 2 - T^Wi = 0, (38) 

2 r 1 ( Di+ ;)*' + ^i L!< ' 2 = a (40) 

The equations above will become more tractable if the variables are changed to 



Then we have 



$o = <t) e lk \ $i = V2r<f> ie lkt , $ 2 = 2r 2 2 e ifc *. 

D + ^ $i -Lx$ =0, (41) 

D - ^ $2 - L $i = 0, (42) 

r 2 B (r) ^D + |^ + M $ + = 0, (43) 

r 2 B(r) (D + -) $! + Lj$2=0. (44) 
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The commutativity of the operators L and D enables us to eliminate each $^ from above equations, and hence we 
have 






L L\+r 2 B(r)[nl + -j I D - ^ 



$o (r,0)=O, 
<I> 2 (r,0)=O, 

$1 M) =0. 



(45) 
(46) 
(47) 



The variables r and can be separated by assuming a separable solution in the form of 

$o (r, 9) = f (r) 9 (9) , $i (r, 9) = fi (r) 9i (9) , $ 2 (r, 9) = f 2 (r) 6 2 (0) . 

The separation constants for Eq. (45) and Eq. (46) are the same, because L„ = — Ljj (n — 6) , or, in other words, the 
operator LjLi acting on 0o (9) is the same as the operator LoL{ acting on 2 (9) if we replace 9 by 7r — 9. However, 
for Eq. (47) we will assume another separation constant. Furthermore, by defining i?o (r) = i ^i( r ) = an( i 



/2W 



the radial equations can be written as 



/o>) + ?/oW+ 



B (r) 



rB (r) B 2 (r) / S 2 (r) r 2 B (r) 



/o(r) = 0, 



(48) 



iuj 



B'(r)\ + ^ 



rB (r) B 2 (r) / B 2 (r) r 2 B (r) 



/ 2 (r)=0, 



(49) 



S(r) 



B 2 ( r ) r 2 B ( r ) 



A(r) = 0, 



where e and 77 are the separability constants. 



1. The case r— > 00 



(50) 



For the case r — » 00, the corresponding metric is given in Eq.(25). Hence, the radial parts of the Maxwell equations, 
(48) , (49) and (50), become 



/>) + ^/>) - 0, i -0,1 



Thus, the solutions in the asymptotic case are 



Co 



Rj(r) = Ci + -£ j = 0,1 
Mr) - C 3 + §, 



(51) 
(52) 

(53) 
(54) 



in which Cj are integration constants. The solution above is square integrable if C\ = C3 = 0. Therefore, the 
asymptotic form of the solutions behaves as Rj(r) <~ ^k, j = 0, 1 and -R 2 (r) ~ 
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2. The case r->- 

The metric near r — > is given in Eq.(27). Hence, the radial parts of the Maxwell equations (48) , (49) and (50) 
for this case are given by 

i<(r) - -i?;.(r) - ^Rj{r) = , j = 1,2 (55) 



Q 
rf 

r " ' ' q 2 



i2o(r) + ^(r)- J^oW = (56) 



whose solutions are obtained as, 

flj-(r) = C 3 e? r (ar-l) + C 4 e-? r (ar + l), j = 1,2, (57) 
R (r) = Q sinh + 9l cosh (58) 

where Cj are constants. The above solution is checked for square integrability. Calculations have revealed that 

Jo 9 2 

which indicates that the obtained solutions are square integrable. The definition of the quantum singularity for 
Maxwell fields will be the same as for the Klein— Gordon fields. Here, since we have three equations governing the 
dynamics of the photon waves, the unique self-adjoint extension condition on the spatial part of the Maxwell operator 
should be examined for each of the three equations. As a result, the occurrence of the naked singularity in f(R) 
gravity is quantum mechanically singular if it is probed with photon waves. 

C. Dirac Fields 

The Newman— Penrose formalism will also be used here to find the massless Dirac fields (fermions) propagating in 
the space of /(-R)-gravity. The Chandrasekhar-Dirac (CD) equations in the Newman— Penrose formalism are given 

by 

{D + e- p)F 1 + (S + TT-a)F 2 = 0, (59) 

{A + n- 1 )F 2 + (5 + f3-T)F 1 = 0, 

(D + e-p)G 2 ~(5 + ir-a)G 1 = 0, 

(A + £-7)Gi- (8 + P-t)G 2 = 0, 

where F\,F 2 , G\ and G 2 are the components of the wave function, e, p, it, a, /i, 7, /3 and r are the spin coefficients to 
be found. The non-zero spin coefficients are given in Eq.(35). The directional derivatives in the CD equations are the 
same as in the Maxwell equations. Substituting non-zero spin coefficients and the definitions of the operators given 
in Eq.(34) into the CD equations leads to 



B(r) 



Do + i) Fl + ^=UF 2 = 0, 



2 \ 2B(r) r) r ^2 
1 




(60) 
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For the solution of the CD equations, we assume a separable solution in the form of 

Fi = /i(r)y 1 (0)e < ( fct+m *'>, (61) 
F 2 = f 2 (r)Y 2 (6)e^ kt+m ^\ 
Gi = gi(r)Y 3 (9)e^ kt + m ^ , 
G 2 = g 2 {r)Y49)e l< - kt+mip) , 

where m is the azimuthal quantum number and k is the frequency of the Dirac fields, which is assumed to be positive 
and real .Since {/i, f 2 , gi, g 2 } and {Y 1 ,Y 2 ,Y 3 ,Y4} are functions of r and 8, respectively, by substituting Eq.(61) into 
Eq.(60) and applying the assumptions given by 

/i(r) = 52 (r) and f 2 (r) = 9l (r) , (62) 
Yl(6) = Y 3 {6) and Y 2 (6) = Y 4 {0), (63) 

the Dirac equations transform into Eq.(64). In order to solve the radial equations , the separation constant A should 
be defined. This is achieved by using the angular equations. In fact, it is already known from the literature that 
the separation constant can be expressed in terms of the spin-weighted spheroidal harmonics. The radial parts of the 
Dirac equations become 



We further assume that 



then Eq.(64) transforms into, 



Do + £)/i(r) = -A=/ 2 (r), (64) 



/i(r) 
Mr) 



fri(r) 
r 

* 2 (r) 



D *! = -^=* 2 , (65) 



r 



V2 



B{r) B' (r) \ A 

2 [ Ua+ 2B(r)J * 2 rV2 V 

Note that \J ^p-T)\\J = Dj + 2B{r) + r ' anc ^ usm S tms together with the new functions 

Ri{r) = *i(r), 



and defining the tortoise coordinate r* as 



Eqs.(65) become 



/?2(r) = tPP*2(r) : 



A +ifc W = ^i? 2 , (67) 

GST"* y 7* 

ik\ R 2 = 

ar* / r 



11 



In order to write Eq.(67) in a more compact form, we combine the solutions in the following way: 

Z+ = R± + i?2 j 
Z— = i?2 — Rl ■ 

After doing some calculations we end up with a pair of one-dimensional Schrodinger-like wave equations with effective 
potentials, 



( — + k 2 )Z± = V±Z±, 



\dr 2 



V± 



BX 



± A— 



d ( VB N 



In analogy with equation (16), the radial operator A for the Dirac equations can be written as, 



(68) 
(69) 



If we write the above operator in terms of the usual coordinates r by using Eq.(66), we have 



d 2 _ B_d_ _L_ 
~ ~~d~^ ~ ~B~dr + B 2 



r z dr \ r 



(70) 



Our aim now is to show whether this radial part of the Dirac operator is essentially self-adjoint or not. This will be 
achieved by considering Eq.(20) and counting the number of solutions that do not belong to Hilbert space. Hence, 
Eq.(20) becomes 



d 2 B^d_ _L_ 
dr 2 + ~B~db- ~ B 2 



=F i ip{r) = 0. 



For the asymptotic case, r — »• oo , the above equation transforms to 

dr 2 r dr 



whose solution is 



ip (r) = Ci + 



C 2 



(71) 



(72) 



(73) 



Clearly the solution is square integrable if C\ — 0. Hence, the solution is asmptotically well behaved. Near r — > , 
Eq.(71) becomes 



whose solution is given by 



(fib 2 dip u 

-tt r- + = 0, 



dr 2 



dr 

a = =F2A<7, 



(74) 



^ W = ( % ) {C3J3 (x) + C 4 N 3 (x)} . 



(75) 



where J3 (x) and N 3 (x) are Bessel functions of the first and second kind, and x = 2y^. As r — > 0, we have x — > 00. 
The behavior of the Bessel functions for real v > as x — > 00 is given by 



(76) 
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thus the Bessel functions asymptotically behave as J3 (x) ~ \j cos ( x ~ t) an< ^ ^ 3 ( x ) ~ y 7rx sm ( x ~ "T") ■ 
Checking for the square integrability has revealed that both solutions are square integrable. Hence, the radial operator 
of the Dirac field fails to satisfy a unique self-adjoint extension condition. As a result, the occurrence of the timelike 
naked singularity in the context of f(R) gravity remains singular from the quantum mechanical point of view if it is 
probed with fermions. 



IV. CONCLUSION 



In this paper, the formation of the naked singularity in the context of a model of f(R) gravity is investigated within 
the framework of quantum mechanics, by probing the singularity with the quantum fields obeying the Klein— Gordon, 
Maxwell and Dirac equations. We have investigated the essential self-adjointness of the spatial part of the wave 
operator A in the natural Hilbert space of quantum mechanics which is a linear function space with square integrability. 
Our analysis has shown that the timelike naked curvature singularity remains quantum mechanically singular against 
the propagation of the aforementioned quantum fields. Another notable outcome of our analysis is that the spin of 
the fields is not effective in healing of the naked singularity for the considered model of the f(R) gravity spacetime. 

Another alternative function space for analyzing the singularity in this context is to use the Sobelov space instead 
of the natural Hilbert space [l3| . The Analysis in Sobelov space entails square integrability both of the wave function 
and its derivative. Although the details are not given in this study, the analysis using the Sobelov space has revealed 
that irrespective of the spin structure of the fields used to probe the singularity, the model considered of f(R) gravity 
spacetime remains quantum mechanically singular. 

Hence, the generic conclusion that has emerged from our analysis is that in the model considered of f(R) gravity, 
the formation of a timelike naked singularity is quantum mechanically singular. 

It will be interesting for future research to extend the quantum singularity analysis in other ETG models. Further- 
more, it will be a great achievement if the criterionproposed by HM is extended to stationary metrics. Although the 
preliminary work in this direction is considered in [3l|, the formulation has not been fully completed. 
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